Constructing explicit magnetic analogies for the dynamics of glass forming liquids 
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By defining a spatially varying replica overlap parameter for a supercooled liquid referenced to 
an ensemble of fiducial liquid state configurations we explicitly construct a constrained replica free 
energy functional that maps directly onto an Ising Hamiltonian with both random fields and random 
interactions whose statistics depend on liquid structure. Renormalization group results for random 
magnets when combined with these statistics for the Lennard- Jones glass suggest that discontinuous 
replica symmetry breaking would occur if a liquid with short range interactions could be equilibrated 
at a sufficiently low temperature where its mean field configurational entropy would vanish, even 
though the system strictly retains a finite configurational entropy. 



The commonalities and contrasts between the glassy 
behavior of supercooled liquids and quenched disordered 
magnetic systems have been long studied. Both exhibit a 
diversity of long lived states with no apparent structural 
long range order. These states challenge the paradigms 
of standard many body physics. Starting in the 1980's, 
the disordered magnetic systems came under control, but 
only in mean field, through the concept of broken ergod- 
icity using replica methods which allowed averaging over 
quenched disorder. Despite the absence of quenched dis- 
order there is an analogy between the first order break- 
ing of ergodicity predicted by mode coupling theory for 
structural glasses and the phase transition predicted for 
spin glasses that lack up-down symmetry (Potts spin 
glasses) [TJ [5] . In mean field, the asymmetric spin glasses 
exhibit a Kauzmann entropy crisis [3] like that for super- 
cooled liquids. 

The analogy between glass forming liquids and mean 
field spin glasses is however incomplete. The observed 
dynamics of supercooled liquids fits the mean field mode 
coupling paradigm only modestly well [4]. The main prob- 
lem with the mean field theory is its neglect of important 
activated motions that allow a supercooled liquid to re- 
configure locally. The most convincing evidence for this 
is the near Arrhenius dynamics with large activation en- 
ergy seen in structural glasses in the aging regime [5J. 

Activated motions have infinite barriers in mean field. 
Describing activated motions requires accounting for the 
finite interaction range. Accompanying a paucity of ex- 
act results even the empirical situation remains con- 
troversial for finite range random magnetic systems jTJ [8]. 
For the random field Ising magnet, arguments based on 
renormalization group theory and droplet arguments do 
work well PI Hul l 11]. but controversy remains as to the ex- 
tent Ising spin glasses partake of mean field versus droplet 
features [IS]. Nevertheless, useful analogies between glass 
forming liquids and disordered magnets have been drawn 
showing how droplet arguments give the Vogcl-Fulchcr 
law [13]. Similar arguments have also been adduced us- 
ing replicas [31 HU [TS] . 



By applying density functional theory to determine 
the parameters in these droplet arguments [16] one pre- 
dicts a large number of confirmed quantitative results 
for liquids [5 . The resulting random first order transition 
(RFOT) theory bears some resemblance to the nucleation 
picture of ordinary first order phase transitions. The the- 
ory's core is an analogy to the random field Ising magnet 
in a field. The average field in the magnet is related to 
the configurational entropy density of the liquid[13 1 (i . 
Tarzia and Moore [17] suggest that supercooled liquids are 
related to Ising spin glasses in a field. When an average 
field is present, corresponding to a finite configurational 
entropy in the liquid, both the random field magnet and 
the spin glass have the same symmetry: neither model is 
expected to have a phase transition although this is still 
somewhat controversial [18 . While both magnetic analo- 
gies agree on this point, the activated dynamics in the 
two different analogical magnets differ because the inter- 
face energies in spin glass droplets have weaker scaling 
than the random field ferromagnet. To quantify this dis- 
tinction in this paper we construct explicitly the analogy 
between a structural glass forming liquid and the corre- 
sponding short range disordered ferromagnet. 

Replica methods along with liquid state theory allow 
an explicit mapping of the free energy landscape of a glass 
forming liquid onto a disordered Ising magnet. Using 
this explicit construction we may eschew droplet argu- 
ments entirely. Simulating the analogous magnet shows 
the mapping reproduces, however, the main features of 
the RFOT analysis. The mean activation barrier and its 
fluctuations which give rise to non-exponential relaxation 
can be computed using importance sampling methods for 
the analog. Approximate arguments can place the analog 
system onto a phase diagram previously deduced for dis- 
ordered magnets using renormalization group (RG) . This 
construction suggests that under thermodynamic condi- 
tions when the mean field estimate for the configurational 
entropy would vanish, the liquid would still undergo a 
phase transition having one step replica symmetry break- 
ing, despite the exact configurational entropy remaining 



2 



non-zero PU [20]. 



I. THEORY 

The analogy yields a description of a structural glass 
in terms of discrete spin-like variables tied to the liquid 
structure equilibrated at one time. A structural glass 
is statistically homogeneous but nonuniform with a den- 
sity /o(x) that is not translationally invariant. Liquid 
state theory provides a free energy as a functional of such 
a density £F[/9(x)]. While the complete equilibrium free 
energy F — — k bT hi Z assumes all phase space can be 
sampled, close to the glass transition there is trapping in 
locally metastable states which manifests itself as exten- 
sively many local minima described by a configurational 
entropy density S c = In N ms . This non-ergodic behavior 
can be captured by a construction due to MonassonpTj 
in which an external random constraining pinning field 
couples to replicas of the system through an attractive 
potential. Each replica's density field is p fc (x). The free 
energy of the m replica system is 



F(m, (3) = lim 
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The typical free energy of a metastable frozen 

v differs from the corn- 
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Monasson's thoroughly field theoretic formulation, 
we separate the to replicas into one fiducial probe 
copy of particles interacting through the Hamiltonian 
H({x{}) = J2i<j u ( x f — x j ) an d to — 1 others described 
by density fields. Here u is the microscopic inter-particle 
potential. The pinning field on the other replicas is 
p^(x) = J^i ^( x — x f)- We can write: 



F(m, 0) = - lim — In / D{xf }e~^ H ^) / J[ Dp fe (x) exp - p £ % fe (x)] 
pm J J fc=1 k=l 
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In calculating the free energy with respect to the probe 
replica we ignore peripheral interactions not involving the 
probe, thus decoupling the partition function 



F(m,0) = - lim J- In / D{xf Ve"^** » [Z q ] m - 1 
5^0 0m J 

(3) 

such that each replica is constrained by a potential to the 
vicinity of the fiducial copy: 



Z q = yDp(x)exp j-/3%(x)] 

■ / dx[p(x)-/(x)] ; 
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For a glassy system, when T5 C = | ro=1 is fi- 

nite, in the limit to — > 1, the replicated free energy is 
dominated by a saddle point corresponding to the spon- 
taneous ordering of replicas in phase space. The saddle 
point solution to the free energy is found by minimizing 
the exponential's argument. 
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The free energy functional can be written, a la density 
functional theory, as the sum of an entropic cost to lo- 
calize the density and an interaction term. 



/3J[p(x)] » J dxp(x) In p(x) + 0Fi, 



t [p( x )]- (6) 



While the free energy functional 3[p\ is globally min- 
imized by a uniform equilibrium solution with mean 
density pa, 7[p\ also has local minima correspond- 
ing to frozen aperiodic densities. The g coupling, 
even as g vanishes, picks out one particular min- 
imum around the structural state p* . Thus ^[p] 
can be analyzed in terms of the similarity, or over- 
lap q = J dx(/?(x) — po)(p * (x) — po), between p(x) and 

p-^(x) = J2i ^( x ~~ x f)- A schematic of the free energy 
as a function of the order parameter q as computed in 
references an d [3] is shown in figure [TJ*. The free 
energy difference between the large overlap solution and 
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FIG. 1: (a) Schematic mean field free energy profiles for 
supercooled liquids at the dynamical crossover temperature 
(dashed line), the Kauzmann temperature (solid line), and 
an intermediate temperature (dot-dashed line). In mean field 
the particle localization, a, and the structural overlap, q, are 
equivalent reaction coordinates. The secondary free energy 
minimum at Tk < T < T c demonstrates the existence of 
metastable structural states in supercooled liquids, (b) Free 
energy profiles calculated for the finite range Ising magnet 
analogous to the LJ liquid. The minimum size needed to es- 
cape the free energy minimum and thus reconfigure the liquid 
at s c = 1.10 is N* = 148 particles. 



the small overlap solution is the excess free energy of 
the frozen glass over the equilibrium free energy and is 
determined by TS C . 

The large overlap state is well approximated by 
a density distribution of a sum of Gaussians cen- 
tered around the particle locations of the fiducial 
(probe) copy p(x) - £,/*(*) = £\ (ff 2 e -^(*-f) 2 . 
The localization parameters, {a/\, determine the lo- 
cal overlap. Near the large overlap minimum, 

q(ai > pl /3 ) = J2i (( Q; i/ 7r ) 3/2 - Po)- In tne opposite 
limit, near the global free energy minimum, the density 



ansatz reduces to the mean density and q({a{} — > 0) = 0. 

For large values of {a;} the particles are local- 
ized very near the fiducial locations {x{} and F 
can be evaluated by using the independent oscilla- 
tor approximation [231 124] which decouples the par- 
ticles at the individual site level. Within this ap- 
proximation F int can be expressed as a sum of effec- 
tive potentials between the interacting density clouds, 
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The localization parameters corresponding to the large 
overlap solution, {aj}, can be found by applying a self 
consistency condition [23J [23]. The existence of the free 
energy minimum at large overlap reflects the cage ef- 
fect where the motion of a particle is restricted by its 
neighbors. To first order one may compute the poten- 
tial PV C [Xi - xf , {aj}) = J2j PVeff (|xj - xjj\;ctj \ and 
expand around small displacements of particle i, Wi — 
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The linear term gives no contribution because the fiducial 
replica becomes centered on a stationary location with all 
forces canceling. The localization parameter of particle i 
thus can be computed from the curvature of the effective 



^1=0, {at}), 



giving a 



potential, aj = gV 2 /3V e (\xi 

self-consistent solution for aj 24, 25J. 

Near the uniformly low overlap state the interaction 
free energy F int follows from the equilibrium liquid equa- 
tion of state, Zeos(v)i where r\ is the packing fraction [25]: 



F^ q = 7Vlnp A 3 - N + N J\z EoS - 1)^ 



(10) 



To characterize the reconfiguration events and de- 
velop the magnetic analogy one must also examine non- 
uniform solutions. At the interface between the two so- 
lutions there must be some energetic penalty due to the 
patching together of distinct configurational states. At 
an interface, one particle's replica is in the large over- 
lap state while it's neighbor has small overlap, so the 

pair interaction becomes flV%^ \\x{ — xj |; a\, at) 

lr J dx i dx J pj(x i )pj(x J )e~5 ,3u ( Xl ~ Xj '). The small over- 
lap parameters {a\} determine p\ and are obtained in the 
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self-consistent phonon theory by matching the entropy 
of the low overlap state calculated within the Gaussian 
density ansatz with the entropic term of the equilibrated 
liquid: 

V^ln^ = iVlnpoA 3 . (11) 
* — ' 2 ire 

i 

For any combination of the discrete values of {a\} and 
{aj}, the free energy of the supercooled liquid is equiv- 
alent to a pairwise interacting model with spins located 
at the fiducial locations, {x{}. 

m=-12 - s *) + 12 J v M 1 - s j) + s j( 1 ^ s ^ > 

i i<j 

(12) 

where the spin, Sj = 1, corresponds with large overlap 
and Si — a small overlap site. The average field is 
found from the bulk free energy difference between the 
states, J^i hi = P F giass - f3Fu q = Ns c /k B , with a het- 
erogeneous local configurational entropy resulting from 
the alpha variations. 

T 

hi = fin ^ +Pl2VeffQx( ~ xj |; a}) - (13) 

The interactions defined through the effective potential 
give the surface energies of droplets within the RFOT 
picture and are explicitly 

Jij = V 2 eff (|xf -xt\;ala})+V 2 ef f (|xf - x/|;aj,aj) . 

(14) 

II. APPLICATION TO A SIMULATED GLASS 

The mapping to the disordered Ising model should be 
carried out for each fiducial equilibrium liquid structure. 
We sample fiducial structures of the Kob- Andersen 80- 
20 mixture of two types of Lennard- Jones (LJ) particles 
at density 1.2 (in LJ reduced units). The pairs have 
interaction parameters [26], [28] . a a A — 1-0, <tbb = 0.8, 
Qab = 0.88, caa — 1-0, £bb — 1-5, cab = 0.5. The fidu- 
cial structures were obtained by simulated annealing runs 
to the temperature Tmd = 0.45. The equilibration time 
for these simulations reaches a tenth of a microsecond 
when referenced to argon. 

The parameters {a}} ({qj}) and {a\} ({qj }) are calcu- 
lated for every particle. The mean RMS deviation deter- 
mined from {cti} is about 0.12 particle spacings, rather 
close to the Lindemann parameter expected for periodic 
crystals g?l ~ 0.1 independent of the force law. The RMS 
actually observed during the MD run at Tmd = 0.45 is 
consistent with this estimate dh — 0.113. 
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FIG. 2: The distributions of interactions and local fields of the 
magnet analogous to the simulated LJ two compound glass. 
In this mapping h is directly related to the configurational 
entropy, h — s c /kB- The fields are shown at h — 1.2, close to 
the dynamical crossover temperature. 



The distributions of calculated interactions are shown 
in figure [2^ as the free energy per neighbor, Ji = 
j: Jij where Zi is the number of neighbors of par- 
ticle %. The typical interaction J = Ji is di- 
rectly related to a, the mismatch free energy penalty 
in RFOT theory for a particle at a flat interface be- 
tween regions of high and low overlap. In the Ising map- 
ping Pa i = ribbJ where ntb — 3.2 is the typical number 
of bonds broken by the interface. The direct calcula- 
tion yields J = 0.55 giving /3a j — 1.77 not very dif- 
ferent from the RFOT theory estimate usually used [2^] 
<Jrfot = ffcsTln^^ = 1.85fc s T. 

In RFOT theory the configurational entropy 
parametrizes a liquid's descent into the glassy 
regime. In harmony with many experimental 
observations |30| l3Tj [32] . the dynamic crossover and 
the laboratory glass transition occur at universal critical 
entropies of s c (T c ) = l.l2k B and s c (T g ) = 0.82k B , 
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FIG. 3: Phase diagram of the Ising model with random bonds and fields adapted from reference [27] ■ The parameters calculated 
for the magnet analogous to the LJ liquid (circular mark) indicate that the liquid would undergo a true phase transition at the 
ideal glass transition. The triangular and square marks indicate estimates of where the glass forming liquids OTP and GeCh 
would fall on the phase diagram. 



respectively. 

Because of the rapidly increasing equilibration time 
scales it is impossible presently to obtain proper fiducial 
structures directly at very low temperatures via molecu- 
lar dynamics. We can, however, treat the configurational 
entropy, and therefore the average field, as variable in or- 
der to extrapolate to find the magnetic system analogous 
to a liquid equilibrated at a much lower temperature, 
eventually extrapolating all the way to the ideal glass 
transition by taking h = -h hi ~ * 0- The presence of 
even a small average field is thought to destroy the phase 
transitions of both spin glasses and the random field Ising 
magnets, but at zero field a transition to a phase with 
long range correlations can still occur. Would a transi- 
tion occur for the liquid analog when h = 0, i.e. when 
the mean field configurational entropy vanishes? We an- 
swer this by appealing to an RG analysis of Migliorini 
and Berker|27 j for the phase diagram for an Ising sys- 
tem in which both the fields and the interactions fluc- 
tuate randomly and independently, encompassing both 
the RFIM and the short range spin glasses. Their model 
is on a cubic lattice. We present their phase diagram 
in figure [3] in terms of the mean field theory based nor- 
malization where both the field fluctuations, Sh, and the 
fluctuations of the interaction strength, SJz 1 / 2 , are nor- 
malized by the total interaction energy per site Jz. This 
parametrization should eliminate trivial near neighbor 
lattice dependence. This zero average field phase dia- 



gram is shown at a temperature T = 1 coinciding with 
the established temperature of the analogous magnet. 
The results for the simulated Kob-Andersen liquid are 
indicated by the dot suggesting the disorder in both the 
fields and the interactions is sufficiently modest so that 
the system would undergo a phase transition to a state 
with infinite correlation lengths and divergent relaxation 
times when the field vanishes. This extrapolation im- 
plies the Kob-Andersen LJ liquid should possess a true 
phase transition to a state with one step replica sym- 
metry breaking (RSB), even though the re-normalized 
configurational entropy, which would include small scale 
droplet excitations, strictly speaking remains finite. For 
such a broken replica symmetry state the free energy 
landscape has divergent barriers between a finite num- 
ber of collective free energy basins each one of which still 
has finite configurational entropy due to local defects. 

According to Landau, the excess heat capacity|33 
yields Ss c — ^/ACpfcg /N corr , where N corr is the vol- 
ume within which the disorder is correlated. The explicit 
mapping for the LJ system gives entropy fluctuations 
5h = Ss c /kB — 0.38, smaller than a typical fragile glass 
forming liquid such as ortho-terphenyl (OTP). Assuming 
a similar distribution of interactions but rescaling the 
field fluctuations to that calculated for ortho-terphenyl 
yields the triangular position on the RG phase diagram 
of figure [3] while the square mark indicates where the 
strong glass forming liquid Ge02 would lie. Even though 
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FIG. 4: (color online) (a) Relaxation times of the Ising model 
analogous to the LJ liquid (circles). The solid line gives relax- 
ation times calculated from free energy barriers. The dashed 
lines show fits using relations derived in RFOT theory (see 
text), (b) The minimum region size able to irreversibly re- 
configure. 



OTP is rather fragile, we see it would still be expected 
to obey one step RSB. 

Our extrapolations neglect any structural changes that 
occur in an actual fluid upon cooling yet we can test 
how well the direct dynamics of the extrapolated analog 
model correspond to droplet analysis. Escape from the 
metastable (s = 1) large overlap state corresponds with 
a large scale activated, structural rearrangement of the 
liquid. Directly simulated escape times are shown as cir- 
cles in figure [4^l. The relaxation time grows rapidly at the 
dynamical crossover temperature appearing to diverge as 
s ci h — > 0. Below the dynamic crossover the growth in 
relaxation time is well fit by lnr/ro ~ s~v . The pro- 
portionality constant for the inverse linear fit (ip — 1) is 
22fce, while droplet arguments in RFOT theory predict 
a slightly larger value In t/t = 32fcs/s c - *p — 2, corre- 
sponding to the unwetted result from RFOT theory [IB]. 



actually gives a closer fit to the relaxation time curve. 
This is consistent with what we have already seen in fig- 
ure [3j that the analog magnet underestimates the disor- 
der in the field moving the system away from the critical 
line where wetting is dominant. 

The average overlap of the liquid frozen density fields 
maps onto the magnetization in the analog magnet, 
q = jfJ2i s i- This coordinate can be used to moni- 
tor escape from local minima. We create free energy 
profiles for this local collective reaction coordinate us- 
ing the weighted histogram analysis method [MJ [35] IMj 
(WHAM). The resulting free energy profile, calculated 
with h = s c /kB = 1-1, is shown as a thick solid line in 
figure [TJd. The metastable minimum at large overlap is 
separated from the global minimum at small overlap by a 
free energy barrier that accounts for the relaxation time 
according to r = r e' 3F * . 

The global overlap is not an ideal reaction coordinate 
for reconfiguration, as it averages over reconfiguration 
events occurring at spatially distinct regions. By select- 
ing a spherical region at random and only permitting mo- 
tion within that region, the overlap becomes a good reac- 
tion coordinate. This is the magnetic analogy of the land- 
scape "library construction" [3T] and is similar to a tech- 
nique recently used to calculate the surface tension near 
a first order transition|37j. By varying the region size, 
the minimum size to irreversibly escape a minimum and 
reconfigure the liquid, N*, can be determined. Free en- 
ergy profiles for several region sizes around N* are shown 
in figure [TJ). The free energy barriers computed for re- 
gions of size N* , converted to relaxation time, are shown 
in figure |4^. Using WHAM and the library construction 
allows descent much further into the glassy regime than 
is possible via direct simulation. The predicted minimum 
reconfiguration size is shown in figure [4]}. At high tem- 
peratures, s c > 1, the growth of region size with decreas- 
ing s c is consistent with N* oc s~ 2 expected from RFOT 
theory, but at low temperatures the growth falls off as the 
result of finite size effects since the cluster size approaches 
the simulated system size itself. The free energy barrier 
is also underestimated for the low temperature range. 

Not all cooperatively rearranging regions are created 
equal. The resulting dynamic heterogeneity of the liq- 
uid is seen in figure [5] showing a collection of free en- 
ergy profiles for different regions at h = s c /kB = 1.1. 
There is clearly a spread of relaxation times which can 
give rise to the stretched exponential relaxation behavior 
4>{t) — e -( t / T ) f>Knn common to glassy systems. If the 
relaxation is entirely heterogeneous the stretching expo- 
nent, Pkww, was shown in reference 129] to be related to 
the spread of free energy barriers, 5F$, through the rela- 
tion (3kww ~ (1 + (8Fi/k B T) 2 y 1/2 . Xia and Wolynes 
argued that regions do not reconfigure completely inde- 
pendently, so free energy barriers larger than the mean 
are lowered by facilitation effects of neighboring regions. 
Using the barrier distribution corrected in this way for 
facilitation yields a non-exponentiality parameter nearly 
independent of temperature, (3kww ~ 0-6, a value char- 
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FIG. 5: (color online) Free energy profiles for different regions 
at h = s c /kB = 1.1. The distribution of free energy barri- 
ers gives rise to the stretched exponential relaxation behavior 
common to glassy systems. 



liquids can be mapped onto a general disordered Ising 
model. This mapping allows a computationally inexpen- 
sive route to low temperature dynamics impossible cur- 
rently by direct simulation. Using this information from 
the replica density functional one can guide molecular dy- 
namics simulations carried out at complete atomic detail 
to more easily reach low temperature structural states. 

Our results suggest the Kob-Andersen liquid should 
demonstrate one step replica symmetry breaking at a 
sufficiently low temperature even though its configura- 
tional entropy including local droplet excitations, strictly 
speaking, will not vanish. The system is, in this sense, 
closer to the random field Ising magnet than it is to the 
Edwards Anderson short range spin glass model. The 
results from the simulation are consistent with droplet 
based predictions using the existing random first order 
transition theory estimates. 
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